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We study here the spherical gravitational collapse assuming initial data to be necessarily smooth, 
as motivated by the requirements based on physical reasonableness. A tangential pressure model is 
constructed and analyzed in order to understand the final fate of collapse explicitly in terms of the 
density and pressure parameters at the initial epoch from which the collapse develops. It is seen 
that both black holes and naked singularities are produced as collapse end states even when the 
initial data is smooth. We show that the outcome is decided entirely in terms of the initial data, as 
given by density, pressure and velocity profiles at the initial epoch, from which the collapse evolves. 
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An important question in black hole physics is whether it would be possible to avoid the naked singularities forming 
as gravitational collapse end states by imposing various possible physically reasonable conditions on the collapsing 
configurations. When a sufficiently massive star exhausts its nuclear fuel, it must collapse endlessly. From such a 
perspective gravitational collapse has been studied within the framework of general relativity extensively for various 
forms of matter such as dust, radiation collapse, perfect fluids, massless scalar fields (see e.g. Ij and references 
therein), to find that both black holes (BH) and naked singularities (NS) result as collapse end states. 

Such collapse studies do assume physically reasonable conditions such as an energy condition, the regularity of 
initial data, and such others. It is however, not ruled out that further conditions, mainly motivated by physical 
considerations, might possibly help ruling out naked singularities, and if this turns out to be the case, this is certainly 
worth exploring in view of the basic importance of this question in the theory and applications of black hole physics. 

Our purpose here is to examine if both BH/NS phases of collapse end states would still occur, when in order to make 
the model physically more relevant, we consider only smooth initial profiles and also include non-zero pressures for the 
collapsing cloud. That is, the initial density and the radial and tangential pressures are assumed to be smooth, C°° 
functions of the radial co-ordinate r. While the Einstein equations by themselves do not impose such conditions on 
the collapse (where the metric and other functions are basically required to be differentiable only) , smoothness of 
initial data is frequently considered to be desirable, especially in view of the numerical evolutions of collapse models 
related to scalar field collapse Also, when the density and pressures are smooth functions, there will be, for 
example, no 'cusps' present initially at the center of the collapsing cloud. Again, while we do not know whether the 
central cores of high density collapsing clouds admit such features or not, it may be a good assumption to make on 
physical grounds, and then to examine the consequences in terms of the collapse end states. 

Our second purpose here is to examine explicitly, how the collapse outcomes are affected when a non-vanishing 
pressure is taken into account. For the case of vanishing pressures, the outcome of dust collapse has been known with 
or without assumptions such as smoothness of initial density distributions |2j ■ Here the initial values of density profile 
and the velocities of the collapsing shells fully determine the collapse end state in terms of either a black hole or a 
naked singularity (see e.g. 4]). The situation is much less explored when pressures are present. We examine below a 
collapse model with a non- vanishing but smooth tangential pressure present. It is then seen explicitly that in fact the 
smoothness of the initial profiles cannot restore the Cosmic Censorship Conjecture (CCC) which forbids occurrence 
of naked singularities in collapse, and that there exist non-zero measure of initial data sets which lead to outcomes 
in terms of either BH or NS. The interesting feature about the model we construct in Section HI is that despite the 
presence of pressure, the final state is determined entirely in terms of the density and pressure gradients specified at 
the initial epoch itself, as in the case of dust collapse. 

The outline of the paper is as follows. In Section II we discuss the collapse equations and the regularity conditions. 
In Section HI, a tangential pressure model is constructed and it is demonstrated how the given sets of initial value 
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parameters such as the initial density and tangential pressure values, decide the collapse final state. The behavior of 
the tangential pressure near singularity is discussed in Section IV, and the nature of singularity is examined in Section 
V. Some conclusions are outlined in Section VI. 

II. EINSTEIN EQUATIONS, REGULARITY AND ENERGY CONDITIONS 

We use the polar coordinates (t, r, 6, (p) to write the spherically symmetric metric as, 

= -e^-^'-^Ut" +e^'^^'^'-Ur'' + R\t,r)dn' (1) 

where d£l^ is the line element on a two-sphere. Also we take the energy-momentum tensor to be diagonal for the 
collapsing Type I matter field (that is, the frame is a comoving coordinate system) which is given by, 

Tl = -p- r; = p.; T^g = = pe (2) 

This is a fairly general class of matter fields, which includes various known physical forms of matter The quantities 
p, Pr and pe are the density, radial pressure and the tangential pressure respectively. We take the matter field to 
satisfy weak energy condition, that is, the energy density as measured by any local observer be non-negative, and so 
for any timelike vector V'^ we have, 

T.kV'V'' > (3) 

This amounts to, 

P>0; P + Pr>0; p + pg>0 (4) 

The dynamical evolution of the system is determined by the Einstein equations and for the metric ^ these are given 
as 

iy'{p + Pr)^2ipe-pr)^-p', (6) 
-2R' + + = (7) 

G-H=l-^ (8) 

where (') and (') represent partial derivative with respect to t and r respectively and, 

Gir,t) = e-^^{R'f, H{r,t) ^ e-^" R^" (9) 

Here F{r, t) is an arbitrary function, and in spherically symmetric spacetimes, it is called the mass function of the 
collapsing cloud, which may be interpreted to give the total mass within the shell of comoving radius r. In order to 
preserve the regularity at the initial epoch, F{ti, 0) = 0, that is, the mass function should vanish at the center of the 
cloud. It can be seen from the equation that the density of the matter blows up when i? = or i?' = 0. Here 
the case R' = corresponds to the shell-crossing singularities. However, it is widely believed (see e.g. 0) that these 
singularities can be possibly removed from the spacetime as they are typically gravitationally weak. Hence, we shall 
consider here only the shell- focusing singularities which occur at i? = 0, where the physical radius of all the matter 
shells goes to a zero value. Now let us use the scaling independence of the coordinate r to write. 



R{t, r) — rv{t, r) 



(10) 
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and we have, 

v{U,r) = 1; f(<,(r),r) =0; v<0 (11) 

where ti and tg stand for the initial and the singular epochs respectively. The condition v < signifies that we are 
dealing with collapse situations only. This means we scale the radial coordinate r in such a way that at the initial 
epoch R = r, and at the singularity, R = 0. The advantage that the introduction of this new variable v offers is that 
the regular center at r = (where we also have i? = 0) is now distinguished from the genuine singularity at i? = 
in that, we now have u — 1 at the initial epoch, and w = at the singular epoch R = 0, but at all other epochs 
in-between v has a non-zero finite value for all values of r. 

From the point of view of dynamic evolution of initial data, at the initial epoch t ^ ti, we now have five arbitrary 
functions of coordinate r as given by, 

v{U,r) ^ va{r)] ilj{ti,r) = i/^Q^r); p{U,r) = po{r) 
Pr{U,r) ^ Pr„{r); pe{U,r) = pg„{r) (12) 

We note that all the initial data represented by the equation H12|l are not mutually independent, as from equation © 
we get, 

Mr)=f{'-^f^-^-^)dr (13) 

Jo \r{po+Pro) PO+ProJ 

Now, to preserve the regularity and smoothness of the initial data let us make some assumptions about the initial 
pressures at the regular center r = 0. Let the gradients of pressures vanish at the center, that is, PrgiO) = Psoi^) ~ 
The difference between radial and tangential pressures at the center should also vanish, i.e. Pro(O) ~P0o(O) — 0- With 
these physical assumptions, and from equation H13|) . it is evident that voir) has the form, 

Mr)=r^g{r) (14) 

where g{r) is at least a function of r at r = 0, and at least a function for r > 0. Hence we see that we have 
a total of five field equations with seven unknowns, p, Pr, pe, ip, v, i?, and F, giving us the freedom of choice of two 
free functions. Selection of these functions, subject to the given initial data and weak energy condition, determines 
the matter distribution and metric of the space-time and thus leads to a particular dynamical collapse evolution of 
the initial data. For collapsing models we have i? < 0. 



III. A COLLAPSE MODEL WITH VANISHING RADIAL PRESSURE AND NON- VANISHING 

TANGENTIAL PRESSURE 

The spherically symmetric collapse models, where the radial pressure is taken to be vanishing, but the tangential 
pressure could be non-zero have been studied in some detail over past years 0- The main motivation here was, 
if pressures are introduced within a collapsing cloud, one would like to understand how the situation differed from 
the dust collapse, especially in terms of formation of black holes and naked singularities as collapse end states. The 
Einstein cluster ^ is an example of such a cloud where tangential stresses are present. This is a spherically symmetric 
cluster of rotating particles where the motion of the particles is sustained by an angular momentum which has an 
average effect of creating a non-zero tangential stress within the cloud. 

What is clear now as we study the models with non- vanishing tangential pressures is that both BH/NS phases do 
develop as collapse end states, that is, the naked singularities of dust collapse do not go away with the introduction 
of tangential pressure. In that sense, these are stable to introduction of pressure. While these are existence results, 
we would like to understand how the nature of the final state will depend on the nature of the initial data given (see 
also 9] for the case of Einstein cluster). 

In the following, we construct an explicit example of a collapse model with a non- vanishing tangential pressure, and 
it is shown that the nature of the singularity in this case can be determined explicitly in terms of the initial values 
of the density and pressure parameters, as given at the initial surface from which the collapse develops. This is then 
parallel to the dust case where such results are available. 

As pointed out in the previous section, we choose the two allowed free functions, F(i, r) and v(t, r) in the following 
way. 



F{t,r) = r^M{r) 



(15) 
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that is, we choose the class of mass functions where is a function of the radial coordinate r only. Also we choose, 

iy{t,r) = c{t) + iyo{R) (16) 

Our purpose here is to examine to what extent the initial data determine collapse final fate, and hence we desire 
to construct here explicitly a specific model of collapse with non-vanishing tangential pressure in which case we will 
show that both BH and NS are possible as end states from smooth initial data, and that it is the initial data set 
which determines the final outcome. This is definitely a generalization over the dust models where initial profiles fully 
determine the collapse end state. Eq. (16), in spite of being a strong assumption, enables us in doing this. 

Now putting equation (|15|l in equation ||SJ), we get, 

3M + r[M,r] 

2? ] 7\~' Pr=0 (17) 

Thus we see that the given choice of the mass function ensures the radial pressure to vanish identically throughout 
the collapse. Also, it is evident that the density at the initial epoch is given by. 

Pair) ^3M{r)+rM{r)^r (18) 

It is clear that, in general, as u ^ 0, p — > cx). Thus the density blows up at the singularity R = which will be a 
curvature singularity as expected. Also using equation (|16|l in equation we have, 

G(t,r) =6(r)e2''"(^) (19) 

Here b{r) is another arbitrary function of r. In correspondence with the dust models, we can write, 

b{r) ^ 1 + r'^boir) (20) 

where bo{r) is the energy distribution function for the collapsing shells. 

Now let us consider a smooth initial data, i.e. the initial density, pressure, and energy distributions are expressed 
as only even powers of 

p{U,r) = Poo + P2r^ + PAT^ + ■■■ (21) 

pg{U,r) ^ + pg^ + ■ ■ ■ (22) 



6o(r) =6oo + &02r' + --- (23) 
With the above form of smooth initial data to evolve the collapse, we can integrate the equation H13|) and get, 

MR) = Pe.R' + ^ . . . (24) 

Using now equation (|16|l in equation ©, we get 

2pg = Ru^Rp (25) 
Finally, using equations H15() . H16|I and (|19|l in equation we have, 

/Ri? = ~a{t)e''°'^^^ V(l + r%o)Re^''° -R + r^M (26) 

Here a(t) is a function of time. By a suitable scaling of the time coordinate, we can always make a{t) = 1. The 
negative sign is due to the fact that i? < 0, which is the collapsing cloud condition. 

Let us define a function h{R) as, 

h{R) = j^, 2g{R) + 0{R') (27) 
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Using equation (|27|l in equation (|26|l . we have after simplification, 



Integrating the above equation, we get, 

t{v,r) 



vv = -y/e'^'^ovbo + e'^"" (y^h{rv) + M) 
y/vdv 



^ye'^'^ovbo + e^^o {v^h{rv) + M) 



(28) 



(29) 



The time of singularity for a shell at a comoving coordinate radius r is the time when the physical radius R{r, t) becomes 
zero. The shells collapse consecutively, that is one after the other to the center as there are no shell-crossings. Taylor 
expanding the above function around r = 0, we get. 



t{v,r) = t{v,0) + r 



dt(v, r) 



dr 



d'^t{v,r) 



r=0 



2! 



r=0 



Let US denote, 



Xniv) 



d"t{v,r) 



dr" 



(30) 



(31) 



r=0 



As we have taken the initial data with only even powers of r, the first derivatives of the functions appearing in above 
equations vanish at r = 0, hence we have. 



Now we can express the next coefficient A2 as, 



Xi{v) = 



1 ^dv[Av^ + Ipey + 2ba2V + f ] 
{hv'^ + boQV + M)i 



where, 



P2P02 



(32) 



(33) 



(34) 



When 1/ = 0, it can be seen from equation (|25|l that the tangential pressure vanishes and the model becomes like 
dust. The time t{v,r) and coefficient X2 for the dust case are given as 



\/vdv 



V y^vbo + {v^h{rv)+M) 



1 Vl^dv[2bo2V + f ] 



V {boov + M)2 



(35) 



(36) 



The important point to note here is that any given initial profile given in terms of the density and pressure values at 
the initial epoch completely determines the function X2 {v) . 
In general, for any constant v surface we have, 



'vv = — vw 



dt 
dr 



We see that the time taken for the central shell to reach the singularity is given as 

^ _ /"^ ^/vdv 

" io y/vbo + v3h{0) + Mo 
From the above equation it is clear that for ts„ to be defined, 

boov + h{0)v^ +Mq>0 
Hence the time taken for other shells to reach the singularity can be given by the expression, 

tsir)^ts, + lr^X2{0) + Oir^) 



(37) 
(38) 

(39) 
(40) 
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IV. BEHAVIOR OF TANGENTIAL PRESSURE NEAR THE SINGULARITY 



From the equations and 1^5)1 , it is clear that near the center the tangential pressure behaves as, 

PS - ^ (41) 

If we wish to calculate the limit of pg along any curve approaching the central singularity then in (u, r) plane along 
all these curve w — > as r ^ 0. Then in {v^r) plane the equation of all these radial curves near the center should 
have the form, 

„ = fcir" (a > 0) (42) 
Rewriting the equation of these curves in (i, r) plane we have, 

dt dt dt f dt\ 



dr dr'dvUr) ^^'^ 

ai Ul UV \«' / along the curve 

Now using (|30|l and calculating the limits on the given curve we finally get in (t, r) plane that the above ingoing curve 
can be written as, 

-t = r2A'2(0) + A;2r(3"/2) (44) 

Now if ^2(0) ^ it can be easily shown that the limit of pe in H41() does not blow up, and the system behaves like 
dust essentially near the central singularity. In case of ^2(0) = it can be shown that for the values of a > 2 the 
pressure blows up. However from the metric it can be calculated that if such curves are to be timelike then strictly 
a < 2. Hence we see that there exist no timelike radial curve along which the tangential pressure blows up at the 
central singularity. However there are a large class of spacelike radial directions, described by ^2(0) = and a > 2, 
along which it blows up at the central singularity. 



Alternatively, we can also look at the solution close to the center. In the limit i? ^ (for the zero initial velocity 
profiles), the equation equation takes the form. 



RR^ = 2 pe^ r^M [r) R^ + r^M (r) (45) 

which can be integrated as 



-^R^F[\^l-^-^R^]=t^U{r) (46) 
^/cj 2 4 4 C2 



where F is hypergeometric function, ci — %pe2 r^M{r) and C2 = M{r). 



Expanding the hypergeometric function F and keeping only the leading order term, we get the solution for R as 



(J^^ rM{r)^{t-ts{r))% (47) 

In the equation l|25(l , if we expand v and p we see that near the central singularity in the limit i? — > the pressure 
pe goes as {pe2Poo){r'^ /R')- Now, in {R,r) plane, along a curve R = kr^ the quantity R' can be calculated from the 
above equation, and when we approach the central singularity we see that the tangential pressure blows up for b > 3. 

We see that close to the center, this curve can be timelike only if 6 < 3. Thus for the curve R — ar^, pg blows up 
only along spacelike directions. A similar analysis was carried out for the curves like t ^ tgo ^ ar^ and it was found 
that pg blows up for the value of f3 for which the curve is spacelike. 
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V. NATURE OF THE SINGULARITY 



We need to determine now when there will be families of future directed outgoing null geodesies coming out of the 
singularity and when there will be none. In the case when such families do exist which terminate in the past at the 
singularity, and which could reach outside observers, then the singularity will be visible. In the case otherwise it is 
hidden within the black hole. Another way to look at this is through the apparent horizon and formation of trapped 
surfaces in the spacetime. As the collapse evolves, if the trapped surfaces form well in advance to the formation of 
the singularity, then the same will be covered. On the other hand, if the trapped surface formation is sufficiently 
delayed during the collapse then the singularity may be naked. The apparent horizon within the collapsing cloud is 
given by the equation , R/ F = 1, which gives the boundary of the trapped surface region of the space-time. If the 
neighborhood of the center gets trapped earlier than the singularity, then it is covered, otherwise it is naked with 
non-spacelike future directed trajectories escaping from it. 

In order to consider the possibility of existence of such families, and to examine the nature of the singularity 
occurring at i? = 0, r = in this model, let us consider the outgoing null geodesic equation which is given by, 

^ = e^-" (48) 
dr 

We now use here a method which is similar to that given in The singularity curve is given by v{ts{r),r) = 0, 
which corresponds to R(ts{r),r) = 0. Therefore, if we have any future directed outgoing null geodesies terminating 
in the past at the singularity, we must have R as t —^ ts along the same. Now writing equation (|48|l in terms of 
variables (u = r"', R), we have, 



du a 



(49) 



Now in order to get tangent to the null geodesic in the (i?, u) plane, we choose a particular value of a such that 
the geodesic equation is expressed only in terms of (■^). A specific value of alpha is to be chosen which enables us 
to calculate the proper limits at the central singularity. For example, for Xi{0) ^ case, we can choose a = | and 
using equation (jHJ, (and considering that i? < 0), we get. 



(50) 

In the tangential pressure collapse model discussed in the previous section we have A'i(O) = 0, and hence we choose 
a = I so that when in limit r — > 0, t ^ tc, we get the value of tangent to null geodesic in the (R, u) plane. 




R> 



dR_3 I R VMoA'2(0) 
du ~7 \ u ) ^/G{VG + Vh) 



(51) 



Now note that for any point with r > on the singularity curve ts(r), we have i? — s- whereas F (interpreted 
as mass of the object within the comoving radius r) tends to a finite positive value once the energy conditions are 
satisfied. Under the situation, the term F/ R diverges in the above equation, and all such points on the singularity 
curve will be covered as there will be no outgoing null geodesies from such points. 

Hence we need to examine the central singularity at r = 0, i? = to determine if it is visible or not. That is, we 
need to determine if there are any solutions existing to the outgoing null geodesies equation, which terminate in the 
past at the singularity and in future go to a faraway observer, and if so under what conditions these exist. Note that if 
any outgoing null geodesies terminate at the singularity in the past then along the same, in the limit as r 0,t tg 
we then have from equation (26) R — 0, therefore H = and G = 1 in this limit as F/R vanishes. Let now xq be the 
tangent to the null geodesies in (_R, u) plane, at the central singularity, then it is given by. 



R dR 
Xq — lim lim — = — — 

t^ts U du 



Using equation (|5T)l . we get. 



(52) 

t^t.;r-tO 



(53) 
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In the (R, u) plane, the null geodesic equation will be, 

R = xqu (54) 

while in the {t, r) plane, the null geodesic equation near the singularity will be, 

t-ts{0) = xori (55) 

It follows that if ^2(0) > 0, then that implies that xq > 0, and we then have radially outgoing null geodesies coming 
out from the singularity, making the central singularity to be a visible one. On the other hand, if ^'2(0) < 0, we will 
have a black hole solution. 

We have, however, already seen in equation that the value of A2(0) entirely depends upon the initial density 
and velocity profiles. Therefore, given any density and tangential pressure profiles of the collapsing matter, we can 
always choose an energy profile so that the end state of the collapse would be either a naked singularity or a black 
hole and vice-versa. 

These features come out clearly in the Figures 1-3, where it is examined how X2 depends upon the various initial 
data values chosen. Here the function ^'2(0) is plotted with different values of initial density, tangential pressure and 
energy profiles. It is seen that there exist distinct regions in the initial data space where a certain region evolves to a 
naked singularity outcome, whereas the other part of the region evolves to a black hole. In the figures the red colored 
surface shows positive values of X2- Whenever the initial data values lie in this region, null geodesies from the center 
come out of the central singularity, as a result creating the presence of a naked singularity as the collapse outcome. 
The green colored region represents black hole. The phase separation is clearly seen. It is to be noted that in case of 
a naked singularity, the singularity curve at the center, as given by equation (|40|l . is an increasing function of r as in 
that case ^2(0) > 0, whereas a black hole solution gives a decreasing or constant curve for shells as the coordinate r 
increases. 

These plots also bring out another important feature arising in collapse, which is that the initial data sets which 
give rise to either a naked singularity or a black hole as final end state are open in the space of initial data from which 
the collapse evolves. 

We note that even though the model considered here has non-zero pressure, this is subject to a strong assumption 
given by Eq. (16). This is equivalent to the ansatz i'{t^r) = h'o(R). A feature that the plots given here bring out 
is, in the space of initial data that is defined under these restrictions the BH/NS outcomes are open sets. Though 
restricted as above, this scenario is more general as compared to dust collapse, in that pressure is included. 

This of course does not mean that, in general, naked singularities would generically develop in any given collapse 
scenario. However, it is to be noted that the issues of genericity and stability are rather involved in Einstein's gravity, 
and are not so well-defined. Hence these must always be studied within the context of the limited models under 
consideration. It is to be hoped that such considerations will possibly throw some light on the nature of cosmic 
censorship statement which one may evolve eventually. 




FIG. 1: Illustration of the dependence of end state of the collapse on initial density po and tangential pressure pe„ for bo2 = 0. 
The red coloured region shows the initial data space which goes to NS, while the green region evolves to BH. 
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FIG. 2: Illustration of the dependence of end state of the collapse on initial tangential pressure psg and energy profile feo2 for 
P2 = 1. The red coloured region shows the initial data space which goes to NS, while the green region evolves to BH. 




FIG. 3: Illustration of the dependence of end state of the collapse on initial density po and energy profile 602 for = 0. The 
red coloured region shows the initial data space which goes to NS, while the green region evolves to BH. 



VI. CONCLUSIONS 



We have considered here gravitational collapse from a smooth initial data. We have also included the presence 
of pressure while developing the collapse models. In this sense, the models here may be considered closer to being 
physically realistic, when compared to, for example, the dust collapse where the pressures are completely neglected. 
Some conclusions and remarks are summarized below. 

1. We see that for the model considered here in Section III, where vanishing radial pressure is assumed and an 
additional assumption (16) is made, both BH/NS phases develop as collapse end states. The important feature that 
we observe is the outcome is fully determined in terms of the initial profiles for density, pressures, and energy functions 
as given at the initial surface from which the collapse develops. This is completely due to the strong assumption 
(16). The point that is made here is, we have a construction or an example to show that naked singularities can still 
develop even when densities and pressures are restricted to be necessarily smooth. Secondly, as seen in the present 
case, the initial data alone could determine the final fate even when non-zero pressures are present in the collaps- 
ing cloud. It remains to be examined if this would hold in still wider classes of models other than those considered here. 

2. We find in the above case, that there is an open initial data space which evolves to a naked singularity, and the 
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same is true for a black hole end state as well. In this sense both these collapse outcomes are generic, within the 

framework of the models considered here. 

3. The current analysis shows that only the requirements of smoothness of initial data, and the inclusion of pressures, 
are not sufficient by themselves to rule out naked singularities. As we see here, smooth initial data does give rise to NS. 

4. Through out we have imposed the smoothness condition, however, it is worth noting that the data is required 
to be smooth only initially. The subsequent dynamical evolution of these initial matter profiles is then fully governed 
by the Einstein equations, and hence at later epochs during the collapse the density and pressure profiles may or 
may not continue to be smooth, and could very well develop non-analytic behaviour. It may be worth examin- 
ing, especially when pressures are included, as to how the smoothness condition may be affected as the collapse evolves. 

5. It is important to note that the analysis presented here only considered locally naked singularities. The null 
geodesies coming out from the singularity can in principle go to infinity, {i.e the singularity can be globally visible), 
or it could also happen that the trajectories could again fall back to the singularity. This in fact depends on the 
global behavior (i.e. for large values of r) of the functions concerned. In the case when a singularity is locally naked, 
one can always choose the global behavior of the rest of the dynamical functions in such a way that the singularity 
becomes globally visible. 
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